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We show through a perturbative and an exact calculation using Wilson's renormalization methods 
that in the problem of two interacting Kondo impurities, on-site potential scattering generates a 
quantum-tunneling between the two impurities through a marginally relevant operator. The magni- 
tude of this tunneling V12 depends on the spin-correlation between the two impurities. For exchange 
interactions between the moments comparable to the Kondo energy Tk, V12 is typically much larger 
than Tk- This implies that the heavy- fermion mass in this interesting range is determined by V12 
rather than Tk- The importance of these results for experiments in coupled Kondo dots is also 
pointed out. 



The two Kondo impurity problem has been of inter- 
est because it presents the competition between mutual 
interactions of local moments and the Kondo renormal- 
ization of the moments to singlet states in a soluble way. 
An aspect of the problem which has received no atten- 
tion hitherto is the the low-energy transfer Hamiltonian 
of quasi-particles between the two impurity sites. In the 
quasi-particle problem for a lattice, this sets the scale for 
the heavy- fermion bandwidth. The issue is also directly 
posed in a recent ingenious realization of the two Kondo 
impurity problem using coupled quantum dots [lj. 

The two Kondo impurity problem was solved by nu- 
merical renormalization group method (NRG) for a 
particle-hole symmetric model obtained by ignoring po- 
tential scattering and was subsequently studied by 
conformal field theory 0, bosonizationQ and other 
methods^. In the single Kondo impurity problem, po- 
tential scattering merely shifts the Kondo resonance with 
respect to the chemical potential in accordance with 
Friedel sum-rule. The sum of the phase shifts in the 
even and the odd parity channels for the two Kondo im- 
purity problem is fixed by the Friedel sum-rule to be tt. 
Without potential scattering the phase shift in both the 
even and the odd parity channels is 7r/2 when the fixed 
point is the Kondo fixed point. Consequently, there is 
no splitting of the even-odd resonances or equivalently, 
no direct hopping of the quasi-particles between the two 
impurity sites. In this case there is a quantum critical 
point (QCP) separating the Kondo singlet regime of the 
impurities from the inter- impurity singlet regime [2j. 

Affleck and Ludwig [3] pointed out that the QCP is 
changed to a crossover by the potential scattering only 
when the exchange scattering at the local moments trans- 
fers conduction electron from one impurity to the other. 
In this paper, we first explain the Affleck-Ludwig result 
by a simple perturbative calculation. We then obtain the 
new result that in such a case an even-odd splitting of the 
resonances is generated which depends on the spin cor- 
relations. We then present [6( exact results of a Wilson's 
NRG calculation to reinforce and quantify these conclu- 
sions by deriving the low energy Hamiltonian. These also 
serve to give the crossover scale of the QCP. The relation 
of the coherence scale to the bare parameters of the prob- 
lem and to the generated spin-correlations is explicitly 



calculated. These results are of fundamental importance 
to the problem of the heavy-fermion state. 

The interactions between the free electrons and the 
local moments in the two Kondo impurity problem are 

-ffimp = 4r (Viofyi + ^2^2) ■ (Si + s 2 ) 
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(Si + S 2 ) 



+ KS 1 -S 2 +J2v(4a^+4A2a) ■ (1) 
cr 

Here Si and S 2 represent two S =1/2 Kondo impurities 
sitting at two lattice sites. "01(2) = Ek^ifzl.ti ^1(2), k are 
annihilation operators for conduction electrons in spher- 
ical waves about the impurity 1 and 2 respectively, a are 
Pauli matrices. J's are the exchange coupling between 
the local moments and the conduction electrons. V gives 
the potential scattering; the maximum magnitude of V 
deduced from the Anderson Hamiltonian is | J|/40- It 
is convenient to choose certain J m so that the generated 
RKKY interaction is while introducing such an inter- 
action K explicitly 0, H[- For the NRG calculation, the 
representation in terms of even (e) and odd (o) parities 
of the conduction electron states is more convenient. The 
latter are obtained by yb e ^ k = (V>i,fc ± '02,fe)/v / 2 while 
Je.o — J+ i J_ . 

Note that in Eq.([T]) the only connection between the 
impurity states is through J_. (The direct hopping of 
the impurity orbitals from one impurity site to the other 
is much much smaller than V, J's or the single-impurity 
Kondo temperature Tk for rare-earth or actinide mo- 
ments.) This alone (i.e. with V = 0) is insufficient for 
quasi-particle hopping between sites 1 and 2. This can 
be motivated through a perturbative calculation by the 
diagrams in Fig. (|T|) . 

The vertex can be in general separated as a vector and 
a scalar part: Ti 2 = Fs,i 2 <5 QQ ' + Tv,i2<?aa' ■ S. Here 
S = Si + S 2 . The scalar part, r,g : i 2 , is an effective po- 
tential scattering electrons between sites 1 and 2. First 
we consider the limit when two impurities are locked into 
either a singlet or a triplet states, S — 0, or 1. The bare 
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FIG. 1: Vertices diagrams for the two Kondo impurity prob- 
lem, (a) shows a bare vertex of J_ while (b) and (c) show 
leading corrections to J_ in particle-particle and particle-hole 
channels, respectively, which contribute a potential scattering 
Vi2iplipi- Solid and dashed lines represent the propagators of 
the conduction electrons and the two impurities with a total 
spin S = Si + S2, respectively. 

vertex [cf. Fig. [T]Ja)] is r[? 12 = J_. This is corrected 
by Fig. QJb) and Fig. QJc) in next leading order. They 
exactly cancel if the intermediate particle-line in (b) and 
hole-line in (c) have the same spectra, i.e. at V = 0, or 
particle-hole symmetry. For V 7^ 0, the generated vertex 
(including similar diagrams with V>ifc as the intermediate 
state) is: 

r s , 12 = P 2 ^±I)j + j_( c y), (2) 

where p is the density of state of conduction electrons at 
the chemical potential, 5 is the total spin of the two im- 
purities and c is a dimensionless constant of 0(1). Such a 
term in two Kondo impurity problem generates the new 
effect, a splitting between even and odd resonances, 

V 12 ^ 2 +h.C. = V 12 (^ e -^ ). (3) 

V12 also depends on the spin-states of the two impuri- 
ties; this is seen from the fact that V\ 2 = V e — V while 
the effective potential scattering including Fig. [Tfb-c) is 

V e ^ e + V ^o 

K,o = V- S{S + 1] cVp\j\ + Jl ± 2J + J_) (4) 

When K -> -00, 5 = 1; then 

Via cx V{ P JJ){pJ + ). (5) 

However, when K — > 00, 5 = 0, the contribution to 
the splitting due to V vanishes. In between these limits, 
there are mixing terms between the triplet and singlet 
states of two impurities, so that the results depend on 
the spin-correlation (Si • Sa). For the general case and 
because at least one of the two vertices in Fig.([T]), J + , has 
singular renormalizations, complete answers can only be 
obtained by exact methods, for example, Wilson's NRG 
method @. 

In NRG, the conduction electron band is discretized 
on a logarithmic grid into a semi-infinite chain with only 
nearest-neighbor hopping terms (here we have two chains 
corresponding to even and odd parity degrees of free- 
dom), while the model can be solved recursively by a 
sequence of Hjy's. 

H N+1 =A 1 / 2 H N + ]T [/ p W P (iv+i) CT +H.c.], (6) 

p=e,o;er 



where A is the discretization energy scale. The impuri- 
ties only interact with the first site (f p oa) through the 
Hamiltonian Hq. From the change of spectra with iter- 
ation, the fixed points and the dimensions of operators 
about them are determined. 

For V = 0, two distinct spectra are identified for large 
iterations with the tuning of K Q • When K is large and 
ferromagnetic, the spectra can be fitted with the effective 
Hamiltonian 

Hn = H N + H ir , (7) 

where the fixed-point Hamiltonians in even/odd iter- 
ations are the same as those of free-electrons in odd/even 
iterations, reflecting the 7r/2 phase shift in both parity 
channels and the formation of Kondo resonances. Both 
resonances develop at the chemical potential, as shown in 
Fig.([2^,). H ir contains terms such as t p fl p fi p +H.c. which 

vanish at N — > 00 as A - ^ -1 )' 2 reflecting the irrelevant 
nature. (t p and other leading irrelevant parameters serve 
to determine Tk)- When K is large and antiferromag- 
netic, the two impurities form a spin-singlet; the eigenval- 
ues of the actual Hamiltonian do not show the even/odd 
interchange, reflecting or 7r (for odd and even parity 
channels, respectively) asymptotic phase shifts. Near 
the QCP at K c w 2.4T£j singular fermi-liquid behavior, 
studied in detail earlier [2, 0] ensues. 




FIG. 2: Spectra for the two Kondo impurity model. The 
location of the resonance with respect to the chemical poten- 
tial is determined from the fixed point Hamiltonian and the 
width from the leading irrelevant operators, (a) shows the 
particle-hole symmetric case; Kondo resonances of both even 
and odd parities appear at the chemical potential, (b) for 
J_ = while finite potential scattering V breaks the particle- 
hole symmetry in each channel; both resonances appear at 
the same energy but tuned away from the chemical poten- 
tial, (c) for both finite V and J_ while K < K c ; the even 
and odd resonances split. For (d), when K is bigger than the 
crossover scale K c , although the splitting increases, the reso- 
nance peaks fall out of the band, reflecting the disappearance 
of the Kondo effect. 

We find in accordance with the perturbative calcula- 
tion above that the fermi-liquid behavior on either side 
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and the singular fermi-liquid behavior near the QCP con- 
tinues for V 7^ 0, J_ =0 but with equal additional phase 
shift in both parities so that the asymptotic eigenvalues 
are given by 



(8) 



with V* — V*. This has a scaling dimension and 
appears as a marginally relevant operator at the two 
strong-coupling fixed points as well as at the QCP. The 
conservation of axial charge 0, H| continues to give re- 
lation between the coefficients for the irrelevant opera- 
tors and they diverge in the same manner as at V — 
near the QCP. When J_ = 0, the two impurities be- 
come effectively two non-interacting Kondo scatters. For 
each channel, a potential scattering induced phase shift 
S v = — tan -1 (irpoV) appears in addition to the ir/2 phase 
shift, as in single Kondo impurity problem [l^]. Here, 
S e = S , two Kondo peaks in local density of states co- 
incide with each other while they are shifted away from 
the chemical potential [see Fig.fl^b)]. 

It is natural to represent the eigenvalues in the even 
and odd parity channels in terms of the phase-shifts S e 
and S in these channels, which can obtained from the 
eigenvalues of the fixed-point Hamiltonian, denoted as 
ui*j for j-th excitation level with p-parity. For free elec- 
trons(with vanishing V and J's), the eigenvalues (77*) can 
be easily obtained numerically for a given A. Their values 
at N — > 00 limit are approximately 

N odd: 7 7 ;=sgn(i)AW-\j = ±l,...,±i(iV + l), 

N even: 77* = sgn^A^- 1 / 2 , j = 0, ±1, . . . , ±^N, (9) 

for 1 < \j\ < N/2. With interactions, the 
eigenvalues^*) can be associated with rj* by a phase- 
shift factor. Quite generally, 



* = A -sgn( I? ; 3 .)'5p/T 77 *. 



(10) 



For instance, at the strong-coupling Kondo fixed point, 
lj* at even iterations (A- 7_1 / 2 A -1 / 2 ) are the same as 77* 
in odd iterations, corresponding to a phase shift S p = 
tt/2. One can then calculate the phase shifts from the 
numerical data following Eq. lfTT)]) . The principal results 
of this paper is the evaluation of S e , S and its dependence 
on KjT^ from which the variation of the coherence scale 
is determined. These results are displayed in Fig. 

In the particle- hole symmetric case (V — 0), S e = 8 = 
tt/2 in the Kondo-resonance fixed point, but S e = n and 
S = in the inter- impurity spin-singlet fixed point. The 
latter merely reflects the fact that in the inter-impurity 
spin-singlet state, the even parity channel is fully occu- 
pied and the odd parity channel is completely empty. An 
abrupt change of the phase shift happens at the critical 
point, K c w 2.4T° , where T° w 1.4 x 1(T 4 L>, is the 
single-impurity Kondo temperature (for J/D = 0.125). 
J_ is then irrelevant. When V 7^ but J_ = 0, S e 
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(a) V=0 

(b) V=0.0125 
J =0 

(c) V=0.0125" 
J =0.005 



FIG. 3: Phase shifts as functions of K. (a)-(c) show the 
cases when V = 0, J- = 0.005; V = 0.0125, J_ = 0; and 
V — 0.0125, J- = 0.005, respectively. J+ is chosen as 0.125. 



and S of two phases are shifted equally from tt/2 and 
7r or 0, respectively. They still behave as two indepen- 
dent Kondo impurities; each gains additional S v from the 
potential scattering term. With both finite V and J_, 
one can observe two features. 1) the phase shift changes 
smoothly at a critical value for K C (V), which smears out 
the critical point. K c is a function of V since the poten- 
tial scattering term renormalizes the Kondo coupling (see 
perturbative calculation above) thus changing the Kondo 
temperature. 2) even and odd phase shifts are different 
while their summation remains 5 e + S a = ir + 2S V , where 
S v can be associated with its counterpart in one Kondo 
impurity problem. In this case, although the axial charge 
in each channel is not separately conserved, their sum- 
mation is. </_ is a relevant operator removing the critical 
point if particle-hole symmetry is absent, V^0. 

Given the phase shifts the potential scattering param- 
eter in the asymptotic Hamiltonian which is a marginal 
operator are determined through S e ^ — — tan _1 (7rpV e * ), 



(V;-V *)(f! feO-flofo ), 



(11) 



which is the coherence operator in effective Hamiltonian 



HI 



Vu^l^+H.c), 



(12) 



where Va = D(V* - V*)(l + A- 1 )/2, i.e, Via in Eq.® 
for A = 1. 

The dependence on the coherence scale V12 on various 
parameters is shown in Figs. (|4l5[) . V12 shows a linear re- 
lation with V and J_ for a whole range of K, as shown 
in Fig.([n^-b). This is in agreement with the perturbation 
results in K — > —00 limit, while as expected its depen- 
dence on J + for a general K is not. The critical point 
becomes a crossover, as the degeneracy of the singlet and 
triplet parts of the Hamiltonian is lifted. Close to the 
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FIG. 4: The dependence of V12 and the Kondo temperature 
Tk (both in unit of the bandwidth W) on the RKKY coupling 
constant K. Other parameters are fixed as J+ = 0.125, J_ = 
0.005 and V = 0.0125. 
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FIG. 5: The dependence of the coherence scale V12 on various 
parameters: a) V; b) J_; c) J+ when K — 0. 

crossover, although V12 increases uniformly, the Kondo 
resonances themselves disappear [cf. Fig.Q]. 

Our most important new physical results come from 
the magnitude of V12 and its dependence on KjT K . In 
the particle-hole symmetric case, such a term vanishes: 
the quasi-particle does not hop directly between the two 



sites. (The analog of this effect in the periodic Ander- 
son lattice model, is that the chemical potential falls 
within the hybridization gap for the particle-hole sym- 
metric case and the system is an insulator 11]). A finite 
V12 has important implications for the heavy-fermion lat- 
tice. Heavy-fermion bands have been derived in a variety 
of ways, all of which have the same physics: constrained 
hybridization of a periodic array of local orbitals with a 
wide conduction band with a matrix element of the or- 
der of the single Kondo impurity resonance width, i.e., 
the Kondo temperature. The heavy-fermion bandwidth 
in such calculations is therefore of O(Tr-). However even 
for Tk — > 0, V12 gives the order of magnitude of the 
splitting of the heavy fermion band at the zone-center 
and the zone-boundary. Therefore, if the lattice prob- 
lem is solved by a self-consistent two Kondo impurity 
problem with a cluster generalization of dynamical mean- 
field theory, the effective bandwidth is then to be deter- 
mined including both the parameter Tk as well as V\2- 
It is therefore important to consider the relative magni- 
tude of these two parameters. The interesting range for 
heavy-fermion compounds is for antiferromagnetic cou- 
pling K/T K of 0(1). In this range for the typical pa- 
rameters shown in Fig. ((4]), V12 is more than an order 
of magnitude larger than Tk- The conclusion then is 
that the effect of V12 dominates in the low energy mass 
and thermodynamic and transport properties. It is inter- 
esting that this term already includes renormalizations 
due to spin-correlations, which are further augmented by 
Landau-parameters . 

In the regime of large K/T K , the Tk given in Fig. (J4j) 
reflects only that the spectral weight at the chemical po- 
tential is only of 0(Tk /V^)- This reflects the disappear- 
ance of the Kondo resonance for such parameters and its 
replacement by the spectra schematically illustrated in 
Fig. (p). 

Recently, double quantum-dot systems have been fab- 
ricated to simulate the two Kondo impurity model [l[. 
When the two dots are allowed to interact, two split 
Kondo resonances appear in the tunneling spectra. The 
spectrum in Fig. (3a) of Ref.[lJ is similar to our Fig. ((5b), 
and the two observed peaks can be identified as the even 
and odd parity resonances of the Kondo dots. 

This work was partially supported by the LANL-UCR 
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